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Abstract 

We consider particles in M. d , d > 2, interacting via attractive pair and re- 
pulsive four-body potentials of the Kac type. Perturbing about mean field 
theory, valid when the interaction range becomes infinite, we prove rigorously 
the existence of a liquid-gas phase transition when the interaction range is 
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finite but long compared to the interparticle spacing. 
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An outstanding problem in equilibrium statistical mechanics is to derive rigorously the 
existence of liquid- vapor phase transition in particles interacting with any kind of reasonable 
potential, say Lennard- Jones or hard core plus attractive square well. This is in marked 
contrast to the situation for lattice systems where proofs of phase transitions abound. Thus 
for an Ising model with ferromagnetic interactions in dimensions d > 2 there are known to be 
two coexisting phases at low temperatures. These are perturbations of the two ground states, 
namely the configuration with all spins up or all spins down. At non-zero temperature, there 
are fluctuations which cause the formation of droplets of the opposite phase, but their energy 
cost is so high that they remain, at low temperatures, in d > 2, only small perturbations of 
the ground state. It was Peierls |l| who first gave a convincing argument of the validity of 
such a picture; the argument was later made fully rigorous by Dobrushin and Griffiths 
@. Independently of this general argument, Onsager Q solved the two dimensional Ising 
model on Z 2 , with nearest neighbor interactions, explicitly and found the behavior of the 
system near the critical temperature which marks the end point of phase coexistence. Since 
that time solutions have been found for many other two dimensional lattice models ||. 
At the same time ferromagnetic and other inequalities as well as the development of the 
powerful Pirogov-Sinai formalism || have resulted in a comprehensive rigorous theory of 
phase transitions in lattice systems, in d > 2, at sufficiently low temperatures. 

The extension of these results to continuum particle systems has proven difficult. The 
ground states of such systems are not at all easy to characterize; they are presumed to 
be periodic or quasi-periodic configurations which depend in some complicated way on the 
interparticles forces. This is however far from proven and hence the analysis of the fluctua- 
tions that appear when we increase the temperature above zero are correspondingly harder, 
indeed very much harder, to study than in the simple lattice systems; moreover key inequali- 
ties like the ferromagnetic ones are no longer available. These problems have been overcome 
so far only for some multicomponent systems with special features. In particular, Ruelle, 
0, proved that the two component Widom-Rowlinson model |5J has a demixing phase tran- 
sition. Ruelle's proof strongly exploits the symmetry between the components present in 
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this model: see also later proofs of phase transitions in related models ||, |10 |, \T%. There 



are also proofs of phase transitions in d — 1 for continuum systems with interactions which 
decay very slowly or not at all. Such models with many particle interactions were analyzed 



by Fisher and Felderhof ||12|| , while Johanson [[13] has considered pair interactions which 
decay as r~ a , 1 < a < 2, proving that at low temperatures the pressure is not differentiable. 

In this letter we report the first proof of a liquid-vapor transition in one-component 
continuum systems with finite range interactions and no symmetries. The basic idea of our 
approach is to study perturbations not of the ground state but of the mean field behavior 
(mfb), i.e. we shall consider situations where the interactions are parametrized by their 
range 7 _1 [|i4| , |T5|| and perturb about 7 = which gives mfb. Such scaling potentials were 
investigated by Kac, Uhlenbeck and Hemmer [KUH] |16[ for a system of one dimensional 
hard rods with an added pair potential 

<t>y(qi,qj) = -"^7exp[-7|gi - qj\), 7,a>0. (1) 

This was later generalized by Lebowitz and Penrose [LP] JDJ to d- dimensional systems with 



suitable short range interactions and general Kac potentials of the form 

7 (&> 1j) = -ai d J(i\qi ~ q-j\) (2) 

with f Rd J(r)dr = 1, J(r) > 

LP showed that in the infinite volume limit followed by the limit 7 — * the Helmholtz 
free energy a takes the form, 

lima(p,7) = CE{a (p) - \ap 2 } (3) 

Here p is the particle density, ao is the free energy density of the reference system, i.e. the 
system with a = in (Q). a is convex in p (by general theorems) and CE{f(x)} is the 
largest convex lower bound of /. (The dependence of ao on the temperature has been 
suppressed.) For a large enough the term in the curly brackets in (^) has a double well 
shape and the CE corresponds to the Gibbs double tangent construction. This is equivalent 



to Maxwell's equal area rule applied to a van der Waals' type equation of state where it 



gives the coexistence of liquid and vapor phases |IT 



Following the work of KUH and LP, various attempts were made to go beyond the 7 = 



limit fTS] , JOJ. It is clear from general arguments, and it follows also explicitly from flU 
that in d = 1 there is no phase transition for 7 > 0. Straightforward expansions in 7 are 
therefore bound to fail for d = 1, in the two phase region. In d > 1 these schemes give 
plausible, but uncontrolled, approximations. The main difficulty comes from the fact that 
the phase transition is a singular event, whose dependence on the parameter 7 is not at 
all smooth. To overcome this problem requires all the modern machinery of Pirogov-Sinai 
theory built up in the past twenty five years [6, 20, 21] plus considerable additional effort. 

It is the success of such an effort which enables us to show for some systems in d > 2, 
that their behavior at finite 7 > are close to mfb at 7 = 0, so that a phase transition 
in the latter yields a phase transition in the former for sufficiently small 7. Such results 
have been recently obtained for Ising models, |[20| , where one uses a version of the Peierls 
argument, exploiting the spin flip symmetry of the model. The absence of symmetries in 
our case requires instead the whole machinery of the Pirogov-Sinai theory ||. To insure 
stabilization against collapse, which would be induced by a Kac attractive pair potential, 
the natural choice made by KUH and LP is to replace point particles by hard spheres or 
similar strongly repulsive pair interactions. Our approach however does not work in such 
a case, as we need a cluster expansion for the unperturbed reference system (i.e. without 
the Kac interaction) at values of the chemical potential or density for which it is not proven 
to hold. Instead we consider point particles and insure stability by introducing a positive 
four body potential of the same range as the attractive two body one. The unperturbed 
system is then the free, ideal gas for which the cluster expansion holds trivially. The price 
is a much more involved mean field analysis, which requires a special choice of the form of 
the interactions. 



We now specify the model, state precisely our results, and give a flavor of the proof |35 
Let q = {qi},i = 1,2,... be a configuration of particles in a domain A C R d ,d > 2. The 



energy of the configuration q is 

#7(?) = -^EE4 2) fe,fe) ( 4 ) 

Here 

4 2) fe,fe)=7^ (2) (79n,7fe), (5) 

■^(fcufta.fcsjfcj =7 M " / " (4) (7?ii ) ---,7&), (6) 

and 

J {2) (r u r 2 ) > 0, J (4) (ri,r 2 ,r 3 ,r 4 ) > 0, (7) 

are fixed, bounded, translation invariant functions of finite range: they vanish whenever any 
of the distances |r» — rj\ is larger than some fixed length l^. 

The equilibrium properties of this system are specified by a grand canonical ensemble 
with reciprocal temperature /3, chemical potential A, and suitable boundary conditions (be), 
i.e. by the Gibbs measure //^ » A . To prove coexistence of liquid and vapor phases for some (3 
and A we have to show that by choosing two different be, one favoring the liquid and another 
the vapor phase, call them + and — , the Gibbs measures obtained in the limit A y M. d , 
describe two phases differing primarily by their densities; the appropriate order parameter 
for this transition. We do this in detail |21| for a particular choice of the interactions 

jW(n,r 2 ) = \B( ri )nB(r 2 )\, (8) 
J (4) (n,r 2 ,r 3 ,r 4 ) = \rf =l B{r J )\ (9) 

where B{r) is the ball in W 1 of volume 1 and center r, i.e. J( 2 \ri, r 2 ) is equal to the overlap 
volume of the two balls (of radius 7r -1 / 2 and (4/37T)" 1 / 3 in d = 2,3) centered at r\ and r 2 . 
Similarly J^'(ri, r 2 ,r%, r 4 ) is equal to the overlap volume of four such balls. 
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Theorem. Let j3 c = and /3 > /3 C f as below), then, for any (3 G (/3 C , /3 ) there exist 
functions 7o(/?) and A(7,/?) such that for < 7 < 7o(/3) the model with and J^> as 
in (|^) has at least two distinct infinite volume Gibhs measures p~p- These measures are 
translation invariant and ergodic (with respect to space translations), with an exponential 
decay of correlations. They have particle densities respectively equal to p 7l( g,- > and 
p 7) p,+ > P-y,f3-- In the limit 7^0, A(7, (3) — > A(/?), p 7 ,/3,± — ► P/3,± and there exist positive 
constants c and 5 such that, \X(j, (3) — A(/3)| + ^ s=± |p 7 ,/3,s — Pp, 3 \ < C7" 5 . 

The reason for the particular choice of the interactions (||) as well as for the appearance 
of j3o in the Theorem are related to the mfb of the system (^) valid when 7 — > (following 
the limit A / IR d ). The mean field equilibrium profiles are functions p*(r) that minimize 
the mean field Gibbs free energy functional, 

FpM = J dr^(\ogp(r)-l) (10) 
- / dr\p(r) - J dr 1 dr 2 J (2) (r 1 ,r2)p(r 1 )p(r 2 ) 

+ i[ / dri " ' rfr 4^ (4) ( r b • • • ,7" 4 )p(ri) • • -pin) 

where p(r) is a test density profile and the first integral gives the entropy contribution to 
the free energy. 

The minimizers p*(r) satisfy the mean field equation bTp^\j Sp(r) = 0. The resulting 
equation is highly non linear and not very much is known about its solutions for general 
and J( 4 ). For the particular choice (§) Eq. (|10"D simplifies to 

*AA(p)= J dr[^-(logp(r) - 1) - XR(r, p) 

~R\,p) + \R\r,p)\ (11) 

where i?(r, p) is the average density over a ball of unit volume in R d centered at r. It is now 
easy to show, using the convexity of the first term, that the minimizers are always spatially 
homogeneous i.e., they corespond to a constant density p > 0. For such a density, the Gibbs 
free energy per unit volume given in (|10|) takes the form 
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f(p) = rVOogp - 1) - Ap - p 2 /2 + p 4 /24. (12) 

The minimizing density will therefore be a solution of the equation 

f\p) =p-Hogp-\-p + p 3 /6 = 0. (13) 

This equation has a unique solution for all A when (3 < (3 C = (3/2) 3 / 2 , while for (3 > (3 C 
there exists a A(/3) such that there are two minimizing solutions pp^ + > pp_ with X(f3) and 
P{3^± the same as in the last statement of the Theorem. In other words, f(p) is convex for 
(3 < (3 C and has a double well shape of equal height for (3 > /3 C , A = A(/5). Moreover, there 
is a (3o > j3 c , given by the smallest value of (3 > (3 C for which f"(pp,±) = 2(j3p / 3 t ±)~ 1 , such 

is positive and dominates the non diagonal ones. 



S 2 T 

that the diagonal part of s rr^k 
When (3 > (3q the second variational derivative of J~p,\<p) is still positive at p = pp^± but the 
diagonal part no longer dominates. The former case is much simpler to analyze and we have 
so far only worked out all the details for that case. 

To prove our Theorem we carry out a controlled Pirogov-Sinai cluster expansion about 
the 7 = mean field state. The first step in this analysis is a "coarse graining" , in which we 
partition space into cubes of size /, with / very large compared to the interparticle spacing 
but small compared to 7 -1 . Given a particle configuration q we call {p x }, x the centers of 
the cubes, the particle densities in each cube. We then show that the measure over the {p x }, 
obtained by integrating out all the other variables is, to within controllable errors, a Gibbs 
measure with an effective Hamiltonian which is essentially a discrete version of the mean 
field free energy functional (0) with p{r) there replaced by p x , and balls replaced by "lattice 
balls" . The important effect of this procedure is that the new effective inverse temperature 
is f3l d . For 7 small enough and / correspondingly large enough, we are now in the right setup 
for the Pirogov-Sinai theory. The remainder terms are exponentially decaying multibody 
interactions. 

The "ground states" of our "lattice system" corresponding to the vapor and liquid states 
are now defined by ensembles of configurations having the p x "close" to the mean field vapor 
and liquid densities pp t+ and pp-- The analysis of this system is conceptually close to the 
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one used in the extension of Pirogov-Sinai theory to continuous (unbounded) spin systems 



developed in [21]. Our analysis is actually simpler than that in |2l|. Instead of using a clus- 
ter expansion which requires dealing with interactions among many Peierl's type contours 
separating "bubbles" of one ground state inside another, brought about by the extended 
range of the potentials, we use a more analytic approach. We show in particular that the 
restricted effective Hamiltonian giving the Gibbs measures of the lattice system correspond- 
ing to the + or — ground state ensembles satisfy the Dobrushin uniqueness condition. We 
then show that this remains true even after the addition of contours to the ground states. 
From this follows the exponential decay of correlations in the liquid and vapor phases, for 
(3 C < (3 < (3q, stated in the Theorem. We expect to prove that similar results will hold even 
at lower temperatures, /3 > (3q, but, as already mentioned, the proof is now more difficult: 
technically, the equation satisfied by the stationary points of the free energy functional (|10D 
is no longer a contraction and the criteria for Dobrushin uniqueness is no longer satisfied by 
the lattice system. 
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